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ABSTRACT 


A satellite for which the orbital rate is commensurate (or rational) 
with the earth’s rotation rate will describe a constant ground-track 
over a whole number of sidereal days. On such a satellite, the small 
orbit- averaged, along-track force from certain longitude harmonics of 
gravity can build up over many commensurate (synodic) periods to pro- 
duce appreciable changes in the semimajor axis of the orbit and longi- 
tude placement of the ground-track. The libratory nature of these 
changes and the specific, harmonics which can cause them for any res- 
onance case are discussed with reference to a circular orbit. General 
formulas are presented for these librations, with typical errors de- 
termined by numerical integration. Actual orbit data from three 24 
hour satellites are used in connection with these formulae to determine 
the resonant earth gravity sectorial harmonics through third order. 
These are found to be 10 6 C 22 = 2.42 ± 0.03, 10 6 S 22 = -1.44 ± 0,03, 
10 6 C 33 = 0.32 ± 0.11, and 10 6 S 33 = 1.18 ± 0.11. It is shown that a num- 
ber of subsynchronous circular resonant orbits offer similar promise 
in easily and uniquely discriminating other low and high-order longitude 
Harmonics, including H 32 , H 43 , H 54 , H 55 , H 76 , and H 9g . 
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THE USE OF RESONANT LIBRATING 
ORBITS IN SATELLITE GEODESY 


by 

C. A. Wagner 

Goddard Space Flight Center 


INTRODUCTION 


Greenwich 


Greenwich 


* Mass 
Anomaly 


Since the beginning of the space age, synchronous (earth-stationary) high altitude satellites 
have been the object of intense investigation because of the unique advantages they offer for world- 
wide communications. More recently, lower altitude satellites with periodically repetitive ground- 
tracks have also been proposed for convenient communications applications. It has been known for 
some time that all satellites with periods that are commensurate with the earth T s rotation period 
can suffer long-term changes in their elements, because of longitudinal variations in the earth’s 
gravity field (Reference 1). This is true since the satellite with a commensurate period and viewed 
in a system of coordinates rotating with the earth, maintains a fixed relationship to the underlying 
longitude gravity field (Figure 1). Resonant orbits as they are used here are synonomous with com- 
mensurate or constant ground track orbits, that 
is, those whose orbital rates are rational with 
respect to the earth’s rotation rate. 

Recent researches into these perturbations 
have been made byBlitzer, Gedeonand Yiounulis 
in the United States (References 2, 3, and 4); 

A. H. Cook and R. R. Allan in England (Refer- 
ences 1 and 5); Morando in France (Reference 
6); Sehnal in Czechoslavakia (Reference 7); and 
Batrakof in Russia (Reference 8). I will treat 
the nature of these orbit changes, in particular 
that of the semimajor axis, from a very ele- 
mentary orbit averaging viewpoint. Restricting 
our analysis to circular orbits permits easy 

derivation of integrable libration equations for the evolution of these orbits under a restricted 
spectrum of longitude gravity harmonics. We will illustrate the power of these amplified effects 
for satellite geodesy by solving (with apparent uniqueness) for two low order sectorial harmonics 
of the earth’s field from observed ground track data on three 24 hour satellites. In addition we 
will show that many other specific harmonics may be as easily discriminated from similar ob- 
servations on other constant ground track orbits. 

Assume that, because of a gravity anomaly, there is net change in an element over the period 
of commensurability. Then, over many commensurate or synodic periods, no matter how small 



Satellite 

[— \ Ascending Equator Crossing — | A.A.E.C. 

(A.E.C.) 

-►Longitude 

r orbits in n 1 sidereal days, to repetition of ground track. 


The vectors above show the magnitude and direction 
of along track forces, F T , due to the mass anomaly. 

Figure 1 —Ground -track of a resonant circular orbit. 
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the net change is, a large change in that element can be generated as long as the geographic 
configuration of the orbit is relatively unchanged. The large, long term change is merely the sum 
of the many small net periodic changes. 

This paper only concerns the long term effects of longitude gravity harmonics on the semi- 
major axis of circular orbits with constant, or nearly constant, ground-track. For a circular 
orbit, without inclination restriction, it is possible to derive a very good approximation of the 
secular change of the semimajor axis by straightforward orbit- averaging of the along-track- 
disturbing forces. The only component of a general disturbing force which can effect a change 
in the semimajor axis of an orbit is the along-track component F x . This is true because the 
instantaneous semimajor axis of an orbit is a function only of its total energy (kinetic and 
two- body gravitational potential). The only force which does work on, and thus effects the total 
energy of the satellite, is the along-track force. 

On a constant ground-track, or resonant orbit, the along-track forces on the satellite due to a 
mass anomaly are seen to vary periodically, in step with the predominant orbit motion (Figure 1). 
Part of this action is analogous to familiar mechanical resonance and can be imagined to lead to a 
buildup in the eccentricity of the orbit (for example, over a long period of time). But unlike the 
usual mechanical resonance case, the orbit average of the along-track forces because of the anomaly 
will not in general be zero. If the orbit-averaged along-track force is not zero, then the orbit 
energy, or semimajor axis will change over the resonant period. A change in the mean motion 
implies a shift or drift rate in the geographic longitude position of the track as the exact commen- 
surability with the earth’s rotation is destroyed. However, a feedback mechanism operates here 
with respect to the longitude position of the orbit. It is evident that the orbit-averaged force will 
depend strongly on the longitudinal location of the track with respect to the disturbing gravity 
anomalies. Over many orbits, the track may shift to a location where the averaged force reverses 
direction; then the action will be reversed. The ground-track may eventually come to rest and be 
driven in the opposite direction. 

At the same time as the ground-track goes through its long-term oscillation, the orbit oscil- 
lates in unison, as the semimajor axis slowly grows and shrinks. In fact, both of these special 
actions in response to the longitude -dependent gravity forces of the earth are parts of the same 
phenomenon, called resonant libration. The proper analogy of this action is the circular pendulum. 
For the restricted case of circular orbits and with the anomalies defined in terms of an infinite 
series of spherical harmonics, it will be shown easily that the defining equation for the evolution of 
the nearly-constant ground-track orbit due to each individual relevant harmonic is just that of a 
circular pendulum. 

As a consequence of the required resonance condition in terms of the longitude harmonics of 
gravity, it will be seen that only a very few of these harmonics can cause this long-term libration 
on a given constant ground-track circular orbit. This fact will prove to be vital to the special 
usefulness of these resonant orbits as geodetic tools. 
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RESONANT ORBIT EVOLUTION 


In order to keep the analysis here as simple as possible without losing the essence of the prob- 
lem, we will temporarily ignore the oblate character of the earth. So long as the orbit is reason- 
ably circular, oblateness merely effects the semimajor axis necessary for commensurability. The 
assumption is made here that the mean motion is entirely due to a spherical earth. The evolution 
of the orbit under the continually acting residual along-track forces is easily calculated through 
the consequent orbit energy and period changes. 

In the absence of perturbations, an orbit has a repeating ground-track if, and only if, the ratio 
of the satellite's mean motion a> to the earth's rotation rate u is a rational fraction. This 

s e 

fraction can be expressed as 


where r is the whole number of orbits the satellite makes in n' sidereal days, the synodic or 
ground-track repetition period. 


The two-body energy E of a satellite and the period T of its orbit are given by 


E = 


(2a) 


and 


T = 


277a 


3/2 


, 1/2 


(2b) 


where a is the semimajor axis of the orbit and y is the earth's Gaussian gravitational constant. 
The first order differences of E and T with respect to the semimajor axis are 


AE 


/x Aa 

H 7 


(3a) 


and 


377a 1/2 Aa 


(3b) 


The work AE done on the circular orbit satellite by the along-track force, averaged over the 
repetition period n', is 


where 


AE = 277- ar F T , 


(4a) 


Ft = -A FTdt ’ (4b) 

n Jo 

and f t is the along-track component of the perturbation force F. Equating the expressions for AE 
in Equations 3a and 4a yields Aa, the change in the semimajor axis each r orbits due to the orbit 
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averaged force: 


Att a 3 r F. 


If the expression Aa (Equation 5) is substituted into Equation 3b, the change of period each r orbits 
due to f t , at, is 


\2n 2 a 7/2 rF T 

AT " ^V2 ' (6) 

It can be easily shown (See Reference 9 for example) that the drift rate k of the ground track 
due to a small period change AT, from commensurate conditions, after n' days, is 


k 


- 277 n ' 


AT 

T 


radians/n' 


days . 


( 7 ) 


From Equations 2b, 6, and 7, the effect of F x is to cause a drift rate after n' days and r orbits 
from commensurate conditions, of 

-1277 2 r n ' F t 

A. = 2 radians/n' days . (gj 

Since this initial orbit change results in only a small deviation from the original ground-track, 
the disturbing forces will be virtually the same over the next r orbits. Thus, the drift rate under 
this disturbance increases linearly from the initial rate by the addition of an increment of the drift 
rate (Equation 8) each r orbits or n' days. The acceleration of the ground-track is thus the same 
as the drift rate at n' days from exact resonance divided by n' days, or as shown in Equation 9a. 


.. -1277 2 m' F t 

K = ^ y ’ radians/(n' days) 2 

= - 12 tt 2 (r/n')F T # , radians/sid. day 2 , (9a) 

where the primed disturbing force (actually an acceleration or force-per-unit-mass) is nondimen- 
sional and is used to simplify the result, and is given by: 


F # 


F 

(fl/a 2 ) 


(9b) 


The coupled orbit expansion rate a under F T is just the semimajor axis change (Equation 5) re- 
written as a differential equation 

a = 4tt a (r/n' ) F t ' , (10) 


where a is in a units per sidereal day. 
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THE DISTURBING POTENTIAL 


The derivation of the orbit evolution equations has been elementary. Equation 10 is merely 
the orbit averaged statement of the classic Lagrange equation for the variation of the semimajor 
axis. The difficult task is to solve for the orbit averaged perturbation force. This is primarily a 
function of the inclination of the orbit and the longitudinal orientation of the ground-track. 

To solve this problem mathematically it is first necessary to obtain an adequate representa- 
tion of the longitude gravity anomalies. For this purpose, it is convenient to use the spherical 
harmonic series representation of the gravity potential. This complete earth potential on the cir- 
cular resonant orbit satellite is expressed as 


t i 

a 


00 n /n 

LL{° 


)" 


P n m ( sin <£) J nm cosm(\-\ nm ) 


(ID 


where 

R 0 = mean equatorial radius of the earth, 
a = semimajor axis of orbit, 

4 > = latitude (geographic), 
x = longitude (geographic), and 
P n m = associated Legendre polynomial (see Appendix). 


The longitude anomalies give rise to those terms where the subscript m, the longitude frequency, 
is not zero. In the following text, we will use the shorthand H nm to refer to the gravity Harmonic 
of order n and longitude frequency m. It is possible to show rigorously that for a circular orbit, 
the orbit- averaged along- track force due to the longitude independent ( m = 0) zonal harmonics, is 


zero (Reference 9). 

To compute the average along-track force 
due to the longitudinal harmonics, it is most 
convenient to write the potential U nm in terms 
of the Keplerian elements (a, I, 8 , n) of the 
circular orbit (Equation 12a; Figure 2). This 
simplifies the evaluation, since the along-track 
direction for this orbit is in the same direction 
as the mean motion, described by the nodal 
argument 0 . Also, use of the Lagrange equations 



Figure 2— Keplerian coordinates of the 
circular resonant orbit. 
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requires that the disturbing potential be in terms of Keplerian elements. Thus 


— j 
a J n 


(■t) 7~! R { F n"P ( I ) ex Pi[(n-2p)0 + m(n-aj e t 

p=0 


(12a) 


where i = f-1, R{ } = real part of { }, m / 0, and 




i n ~ m (n + m)! 
2 n p! (n-p)! 


£Vl) k COS3—2P-2K (I) . sin „-n + 2p + 2K(I) . 

K= 0 


(12b) 


Apparently, this type development of the longitude harmonics was first used by W. M. Kaula in 
1961 (Reference 10). The specific form above, especially the inclination function F( I ) , and the de- 
velopment below, most closely follows the development by R. R. Allan in 1965.* 

The along-track force on the circular orbit satellite is the partial derivative of the potential 
U nm with respect to the nodal argument theta (Reference 9). The simple orbit- averaging of this 
derivative provides all the necessary insight into the requirements for and conditions of librational 
resonance of the circular orbit. From Reference 9, the orbit averaged force F x ' is given in terms 
of the orbit-averaged derivative by 


a^nm 


f t' = u SS 


(13) 


The partial derivative of the potential U nm (Equation 12a) with respect to 0 , in terms of real phased 
harmonics, with 6 and time as the only variables, is 


dU 

nm 

~w 


J nm (~r) 2Z F n' mp (!) (n- 2p)cos [(n- 2p)0 +m(fi-^ e t~\ nm ) + 8„ m ] 

P=0 


(14a) 


where 


s 


run 


7 T 

2 


(n - m + 1 ) 


(14b) 


and 


F nL(I) 


F (I) 

nmp V / 
£ n~m 


(14c) 


* Allan, R. R., "Resonance Effects for Satellites with Nominally Constant Ground-Track,” R. A. E. Internal Document (Private Com- 
munication), Royal Aircraft Establishment; Farnborough, England, August 1965. 


6 


The orbit-averaged derivative, as a series of n + 1 integrals of simple phased harmonics in the 
time and 6 for each longitude harmonic, is 


au 


36 



dt 



(I) • (n - 2p) cos 


j^(n - 2p) 6 + m(n- < 



(15) 


The nodal argument can be converted to time through the mean motion and the argument at epoch 
( 6 o, when t = 0): 


0 - 6q + t 


(16a) 


Let 


C lp - ( n -2p)0 o + 


(16b) 


and 


C 2 - (n-2p) - mu , (16c) 

p 

where the phase constants C, are functions of the longitudinal position of the orbit. Sub- 
stitution of the expressions in Equations 16 into Equation 15 yields 


<9U 

nm 

~3W~ 


c n — 

/R \ n I n 

= I ^v) F ""p < I )( n " 2 P) cos ( c i I , +C 2 P 4 ) 

Lp =0 


dt 


(17) 


The frequency constants C 2 yield the critical resonance requirements. Integrating the n + 1 har- 
monic terms in Equation 17 over n' days, for the orbit-averaged force, gives 


d U 

nm 



( n ‘ 2p)F ' 



+ n'C 2 

p 


) 


sin C 


,] 


(18) 


If the frequencies C 2 are not zero (Equation 19), but are arbitrary, then we will have a zero result 
for each integral over a time interval (n' )* not greater than one day for r > n ; and n'/r days for 
r<n’ . In this case, the long-term effects on the mean motion due to that longitude harmonic term 
can be expected to be negligible, at least over the maximum such averaging interval. 
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The condition for zero orbit-averaged force is 


(n')*C 2 = +277, ±4t r, ••• • ( lg ) 

However, when the frequency constant is zero, (Equation 20), each term of the averaged force is in 
an indeterminate form for any averaging period. 

The condition for resonance is 


( 20 ) 


By I/Hospital's Rule, this form can be evaluated with the nonzero result 

H.: (-4 
n 'Co ~*o \ n C, 
p \ ■ 

The implication of the nonzero averaged force, or resonance, condition, from Equations 20 and 16c 
is: 


;in (c. +n'C 2 \ - sin Cj 

V p p / 


cos C 


( 21 ) 


C 2 p - ^ s (n- 2 p)-m cu e - 0, (22a) 

or 

"• _ m _ r . , 

— - ’ a rational fraction . (22b) 

Hence, as expected, the nonzero averaged force occurs only for commensurate, constant ground- 
track orbits. 

As stated previously, with the rational fraction r/n' reduced to lowest terms, n' will be the 
number of sidereal days over which the ground-track of the satellite will repeat (after r orbits). 

In addition, this single resonance condition (Equation 22) tells what harmonics n, m will be active 
on a given commensurate orbit specified by r andn' (Equation 23). 

On given resonant orbit (r/n'), only those harmonics H nm are active for which 


m - ^7 (n-2p) ; p - 0, 1, 2, • (23) 

Furthermore, for a given resonant orbit r/n' and active harmonics n, m, only one of the n + 1 force 
integrals is relevant (nonzero over the averaging time n'); that integral is the one whose p subscript 
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also satisfies the single resonance criterion (Equation 22), where 


i 


P = V = °> 1. 2 - 


( 24 ) 


For the determined r, n ', n, m, and p of the resonant orbit-harmonic, n' C 2 =0. The single nonzero 

p 

force integral is (from the result of Equation 21 in Equation 18): 


dU~ 

nm 

~dF~ 


^ 

/ R o\ mn' fmn' , . "1 

-r F n m p(I)COs[~ «o+m(n-\n.) +8 0m J 


(25) 


where 


n - 2p - — • 


ORBIT EVOLUTION EQUATIONS 

Let k be the geographical longitude of the ascending equator crossing just prior to the epoch 
(Figure 3). Then, assuming fi and a> s constant, since, for the resonant orbit co e t' = 0 O (n'/r) : 

\ = T + n • (26) 

It is recalled that, in our formulation of the 
potential (Figure 2), epoch is the moment when 
the Greenwich meridian crosses the vernal 
equinox. 

From Equations 13, 25, and 26 we finally 
have the averaged force F T ' for each relevant 
harmonic as a function of the inclination and 
longitudinal position of the ground-track, and 
other variables (Equation 27). The relationship 
between the averaged force and the acceleration of the ground-track having already been estab- 
lished (Equation 9), the desired relationship of this acceleration to the longitudinal position of the 
nearly constant track is given in Equation 28, which is a pendulum equation for each relevant 
harmonic. Thus: 



J nm (ir) “ {IF) Cp C 1 ) COS [ m (^ “\.J + 8 om] 5 

n,m relevant 


(27) 
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* 


12 n 2 


L 


m F ' (I) cos 

nmp v ' 


[ m (\- X nm ) + S nm ] radians/sid. day 2 


n ,m relevant 


(28) 


The coupled long-term semimajor axis variation under f t ' is found by substituting Fy from Equa- 
tion 27 into Equation 10. More generally, if l = 1 refers the track to the longitude of the descend- 
ing equator crossing, and 1 = 2 refers it to the ascending equator crossing, the parallel long-term 
motion of all the equator crossings can be derived from 


H - 1 2 tt 2 2 ^ JnmVr) mF nmp( i )c° s [ m (M ~ +S ' ( n - m - *)] radians/sid. days 2 , (29a) 

n,m relevant 

where 


S' (n, m, %) - S nm + 7r4(n -m) 


= \ [(n -m) (2^ + 1) + l] . 

The coupled semimajor axis variation with respect to either node is then 


(29b) 


~ 4tt a 


( R \ n 

_ -t) mF ™ p ( i ) cos [ m (^~0 +S ' ( n - m > ^>1 (a units )/sidereal day 


(30) 


n,m relevant 


While the longitude regime equations (Equations 29a and 30, see also Reference 2) apply only under 
commensurate, circular-orbit conditions, they have been shown to give a very good approximation to 
the resonant orbit evolution for moderate equator crossing drift rates and orbit eccentricities not 
exceeding about 0.01 (Reference 11, 12). 

In reference to present knowledge of the earth’s gravity field, the compound "pendulum 
drift" equation and its coupled semimajor axis variation appear to describe very closely the 
long-term regime of all near-circular, near- resonant orbit satellites under the action of the 
longitude harmonics of gravity. This is true even for longitude drift rates well into the 
circulating-pendulum or world circulation regimes. The numerical studies in References 9, 

11, and 12 show that longitude accelerations estimated by Equation 29a should be accurate to 
within 1 percent of maximum accelerations for any resonant regime. This accuracy should 
hold so long as the drift rate is not substantially greater than the maximum permitted for 
libration (Table 1). The eccentricity of the orbit should also be less than 0.01 to assure this 
desired accuracy. 
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Table 1 


Libratory Characteristics of Geodetically Promising Resonant Circular Orbits.* 


Orbit period 
(Sid. hours) 

Approximate 

semimajor 

axis 

(earth radii) 

Clearly 

dominant 

resonant 

harmonic 

(H „„> 

Inclination 

(degrees) 

^"max 

(10' 5 rad./day 2 ) 

k 

max 

(deg. /day) 

T . 

nun 

(years) 

T /4 

min' 

(days) 

Maximum 

libration 

excursion 

(360/m)° 

(degrees) 

3.00 

1.653 


60 

0.0752 

0.035 

— 

7.01 

640 

45 

4.00 

2.002 

H ,6 

60 

0.0753 

0.041 

3.11 

740 

60 

4.80 

2.261 

H SS 

60 

0.2462 

0.080 

4.90 

447 

72 

6.00 

2.623 

H 54 

60 

0.1025 

0.058 

8.43 

775 

90 

8.00 

3,178 

H 33 

60 

1.086 

0.22 

3.02 

275 

120 

8.00 

3.178 

h 33 

30 

0.4502 

0.14 

4.69 

428 

120 

12.00 

4.164 

H 3 2 

30 

0.3481 

0.15 

6.52 

595 

180 

14.40 

4.702 

H 5S 

30 

0.0060 

0.0125 

31.4 

2860 

72 

14.40 

4.702 

H 55 

60 

0.0096 

0.0158 

24.8 

2270 

72 

16.00 

5.044 

H 43 

30 

0.0342 

0.038 

17.0 

1550 

120 

24.00* 

6.61 

! 

h 22 

0 

2.92 

0.438 

2.25 

205 

180 


•All libratory data, except that for the 24-hour satellite, is derived from a planning curve derived by W.M. Kaula (Reference 15), (Figure 4). 
*Froin actual data (Reference 11). 


LIBRATORY CHARACTERISTICS OF RESONANT ORBIT SATELLITES 

From Equation 29a, we can see that if the inclination and semimajor axis undergo only small 
changes, the long-term longitudinal drift of a nearly-constant ground-track orbit because of any 
relevant harmonic H nm will be analogous to that of a circular pendulum. In fact, many integrations 
of the complete equations of motion have confirmed this assumption for a wide variety of resonant 
and near- resonant orbits (References 13, 14, 5). As such, the drift as a function of time (or, more 
precisely, the time as a function of drift) can be given exactly only in terms of an elliptic integral 
(and then only for a single harmonic effect), (References 13, 5). However, the drift rate regime and 
the critical or equilibrium points of the pendulum drift can be found directly from Equation 29a. 
Additional information about the libratory drift will be the minimum oscillation period (about 
the stable equilibrium points) and the maximum drift rate possible for "capture" in the libra- 
tory regime due to any relevant harmonic. This information also is easily found from Equa- 
tion 29 a. 

The standard pendulum equation for the drift due to H nm will be taken as: 


+ Kl sinm(\-\ e ) - 0. 


( 31 ) 
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where, for convenience, \ refers to the ascending equator crossing longitude; \ e is a position of 
stable equilibrium (k e s ) if K n 2 m > 0, and unstable equilibrium (\ e _ u ) if K n 2 < 0. These positions are 
separated by 


\ - \ 
e , s e ,u 


Equation 29a can be put in the form of Equation 31 as 

= K nm COS [4'^nJ + m (\, m "O + Ti] 


(32) 


(33) 


where 

K „l = 12772 Jm.^) mF nrnp( i ) rad./sid. day 2 . 

Since J nm < 0, the positions of unstable equilibrium are given from Equations 33 and 29 by 

= S'(n, m, 2), for ^nmp > 0 . 

or, (using Equation 29b with l = 2 and the fact that n - m is an integer or zero), 


(34) 


.. ‘ 2S ( n_m ) 


(35) 


Similarly, those positions of stable equilibrium for F n ' mp > 0 are given from Equations 35 and 32 by 


^e,s - 2E [ 2_ ( n_m )] + K m • (36) 

It may be noted that for low -inclination satellites (i '< 30°), F n | np is always positive, and for these 
orbits, Equations 35 and 36 will be the equilibrium criteria. For example, for the 24-hour 
satellite, the relevant H nm are all those for which n-m is even (see Equation 23), orn-m = 2K; 

K = 0, 1, 2, • • •. Thus, for the low-inclination, 24-hour satellite, the unstable equilibrium ascend- 
ing equator crossing longitudes are at 


,u ^nm m 

and the stable equilibrium longitudes for low-inclination, 24-hour satellites are at 

V, = K. + ©O-K) - 


(37) 


(38a) 


where 


K 


= 0, 1, 2, 


(38b) 
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The situation is just the reverse for inclinations (at medium and high inclinations, for example) 
where F n ^ p <0. In these cases, the stable equilibrium longitudes are given by Equation 35, and the 
unstable by Equation 36. Thus, 


5S ( n_m ) 


(39) 


and 


= A. 


sr [ 2 ~ ( n_m )] 


(40) 


for resonant orbits where F ' < 0. 

nmp 

Since the unstable and stable positions are equally spaced around the globe and separated by 
77 /m radians, there will be a total of 2 m of these around the equator. There will thus be m unstable 
longitudes and m stable longitudes for each relevant H nm . The maximum libratory excursion of the 
ground-track will thus be 360°/m. 

Proceeding with the integration of the pendulum equation (Equation 31), the variables ( A .) 2 and 
X. can be separated by means of the differential reduction 


X = = . 

2 \ dt 2dK 

First Equation 31 is rewritten to apply to \ e = k e s : 

k + Ik * I sinm (\ - \ ) = 0 . 

| nm j \ e ,s J 

Substitution of Equation 41 into Equation 42 separates variables: 


(41) 


(42) 


d(X ) 2 = - 2 |k „ 2 j sinm(X -\ e> .) dX 


(43) 


The integral of Equation 43 is the energy, or first integral, of the pendulum equation: 

2 1 R n 2 | 

(^) 2 3 C o + S - c °sm(X-X e s ) . 

If the initial condition X = X 0 ,atx = X 0 ,is given, then Equation 44 may be written as 


(44) 


<M a = (^o) s 


2 Kj 


[cosm(X-X e s ) -cosm(X 0 -X e j] 


(45a) 
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From Equation 45a can be found the maximum drift rate change possible in a resonant orbit due to 
each relevant harmonic. This will occur with nearly exact resonance (\ 0 = 0) momentarily over 
a position of unstable equilibrium. Then the drift rate passing the stable point ( 77/m radians away 
from k eu ) will be the maximum possible libratory drift rate. For k Q = 0, \ 0 = k e u = \ e>s “ 77 /m , 
and Equation 45a reduces to 


2 Kj 


(M 2 


[cosm(\-\ es ) +l] 


The maximum libratory drift rate (at k = \ ) will then be 


(45b) 


(MmAX,™ 



(46) 


Since the maximum acceleration in the H nm libratory regime, from Equation 42, is |\ MAXnm | = [k^ 
Equation 46 can be written as 




2 



MAX ,nm 


m 


(47) 


Another parameter of interest with respect to the regime of the resonant satellite is the mini- 
mum libration period. This is the period of small angle "pendulum" oscillations of the ground-track 
about the stable equilibrium points defined in Equations 36 and 39. If k -k e is always sufficiently 
small (as in a simple pendulum oscillation), sinm(\ -\ e s ) can be replaced by m(\ -k e s ), and 
Equation 42 becomes a simple harmonic equation 


* + | K n 2 J = 0. (48) 

From Equation 48 we can find the small oscillations or minimum period libration about the stable 
equilibrium point. This minimum period is 

_ 2E _ 277 _ 277 

T MIN,nm Ct> f~\~ T7 " 7 I? T * / Af\\ 

y m | K nm| y m |^n m ,MAx| ^ 

Table 1 gives these three libratory characteristics, X MAX , (\) MAX , andT MIN , for the most geodetically 
promising (underlined) orbits (see the next section). A comparative case, from actual data on 24- 
hour satellites, is provided as the last example in Table 1. The scale of these librations can best 
be appreciated by a study of the data in the last two columns of this table. The full librations of 
the ground-tracks have periods of years or tens of years, and amplitudes of tens of degrees. 

The complete analytic treatment of the drift regime due to a single resonant harmonic (for a 
circular orbit) can be implemented by solving Equation 45 in terms of elliptic integrals. This will 
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not be done here. The form of the general solution has been shown previously, (Reference 13), as 
has the specific solution for the long period of libration of one day resonant orbits.* As in the 
case of the circular pendulum, the actual libration period is strongly dependent on the librational 
amplitude. In addition, Gedeon et al.* has also obtained analytic solutions for the 1-day librational 
resonances of an eccentric orbit due to single harmonics. 

We note only that, since the basic averaged equation of librational motion, Equation 29a, is 
nonlinear, the complete solution of the quasi-libration due to the total, infinite set of resonant har- 
monics, cannot be a simple superposition of the individual harmonic solutions. In fact, there is no 
second integral to Equation 29a which is in a simple closed form. However, the first, or energy, 
integral to Equation 29a can be written down analogously to Equation 45a. In terms of i = 2, or 
the longitudes of the ascending equator crossing, Equation 29a may be rewritten with the differen- 
tial reduction (Equation 41), as 


d (X) 


M 2 = 2dX £ K„’ cos[m(\-\ n J + S'(n,m, 2)] 


n,m relevant 


(50) 


The integral of Equation 50 in terms of the initial conditions k = A. 0 at k - k Q (t = 0) , is 


(M 


2 - 


^ K 2 r 

( ^o) 2 + 2 { sin [ m ( X " ) + S ' ( n - m - 2)] “ sin [ m ( A -o “ \.) + S ' ( n - m - 2)]| • ( 51 ) 


n,m rel . 


Then, formally, variables in Equation 51 can be separated to find the time as a function of drift: 


* s 

dM(\,) 2 + 2 YL Sin (' 

- \ nm ] + 8 * [n, m, 2]^ - sin 

+ s< [n - m - 2] )] 

^j-l/2 

• (52) 

J\ C n.mrel. 

/v n 



J 



In all of this simple derivation, the elements of the orbit have been assumed to be essentially 
constant (in particular K nm ), not only over the resonance period n' , but over the long librational 
time t . While this will undoubtedly be a very good assumption for the semimajor axis (which in 
most cases will suffer small relative changes over long periods of time), it may not be so with re- 
gard to the inclination. Because of the gravitational action of the sun and moon and the earth zonal 
gravity, the orbit plane of a satellite (otherwise unperturbed) will precess about a plane between 
the ecliptic and the equator, causing a long-period sinusoidal-like change in the satellite’s in- 
clination (Reference 16), The period of this inclination change (with an amplitude of the order of 
about 10 degrees) can be as low as 10 years for close-earth satellites. However, for the close- 
earth satellites (a* < 4, for example) zonal gravity equatorial precession predominates, so that the 


*Gedeon, G. S., et. al., "Orbital Resonance Caused by Tesseral Harmonics,” TRW Systems Internal Document (Private Communication); 
Redondo Beach, California, (Sept., 1966). 


15 



total inclination change is small. Thus, it will probably be found that only the relatively weak, far- 
earth resonances can be disturbed seriously by these plane precession effects. 


SENSITIVITY AND DISCRIMINATION OF LONG-TERM HARMONIC EFFECTS 
ON RESONANT ORBIT SATELLITES 


For any given commensurate orbit, only a limited number of harmonics cause these resonant 
librations (Equation 23). On all but 12- and 24-hour satellites having resonant orbits at or below 
synchronous, only two harmonics of low order n are active, and the lowest is separated from the 
next highest by two orders of n . Since the strength of gravity harmonics declines as the inverse 
n th power of the distance from the mass center of the earth, it is reasonable to assume that for 
many of the resonant orbits, especially the more distant ones, only the single lowest-order longi- 
tude harmonic will dominate the libratory regime. 


If this assumption is valid, then only two well-defined and separated measurements of the lon- 
gitude acceleration of the nearly- constant ground-track will be necessary to define the amplitude 
and phase angle of these dominant harmonics. To find those resonant orbits which appear likely to 
show the strongest librations and to have the most favorable damping effect on higher order har- 
monics, a curve for the decline in the strength of the harmonics has been used, which is derived 
from a recent study of satellite data by W. M. Kaula (Reference 15 and Figure 4). This study re- 



Figure 4— J nm as a function of n 
from recent satellite studies. 


ports so called normalized harmonic coeffi- 
cients. The relative physical significance of 
the gravity harmonics may be found by a direct 
comparison of these normalized coefficients. 

It can be shown (Reference 17) that the nor- 
malized harmonics, P n m (sin <t>) cosm(\ -\ nm ), 
all have a mean square amplitude of unity over 
the surface of the geoid (actually a sphere of 
unit radius) if = -1. The relation of nor- 
malized to unnormaiized coefficients is shown in 
Figure 4. The individual gravity coefficients which 
affect the resonant orbits, have uncertainties 
which encompass the curve in Figure 4 in almost 
all cases. This curve apparently reflects the fact 
that the strongest longitude gravity anomalies 
are rather deep in the mantle of the earth. It 
has been accepted provisionally only for plan- 
ning purposes. Using it, a table of dominant 
and subdominant, or next strongest, resonant- 
acceleration effects on a large number of sub- 
synchronous orbits has been calculated from 
the longitude regime equation (Table 2). The 
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Long Term Geographic Accelerations on Longitude Gravity Resonant Orbits. 
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/:Checked effects are reasonably strong, but use of such single orbits probably will not permit unambiguous discrimination of these harmonics. Reasonable strength for these 
effects are based on the tracking record of "Early Bird” (Reference 11). [sensing 0.01 x 10 5 rad. /day 2 in 60 days] 




maximum harmonic accelerations in these orbits occur when the equator crossings of the constant 
ground- track are midway between one of the m stable longitudes and one of the adjacent m unstable 
longitudes due to a given harmonic. The accelerations are strongly dependent on inclination, through 
the sinusoidal inclination functions in the regime Equation 29a. In Table 2, those longitude har- 
monic accelerations which seem to offer the best chance to be quickly and unambiguously dis- 
criminated have been underlined. By best chance, we mean they appear to be resolvable in 
a few month 1 s time, and of an order of magnitude stronger than the next strongest harmonic ef- 
fect. This criterion for quick resolution of these effects is based on the tracking of "Early Bird" 
(Reference 11), which resolved an acceleration of 0.01 x 10 5 rad/day 2 in 60 days. 

For 30-degree inclined satellites, it seems that the resonant orbits of 8, 12, 14.4 and 16 hours 
are especially suited for discrimination of the harmonics H 33 , H 32 , h S5 and H 43 . At higher in- 
clinations and lower altitudes, the higher-order longitude effects are stronger and stand out more 
clearly. For example, orbits of 3, 4, 4.8, and 6 hours with about a 60-degree inclination appear 
favorable for discrimination of the harmonics H 9g , H 76 , H 55 and H S4 . 

The observations necessary to detect these long-term amplified gravity effects are for the most 
part elementary and easily implemented. The orbit at the start of a free-drift period must be reason- 
ably circular and near-stationary. The inclination should also be known to within a few tenth’s of a 
degree. During the long free-drift periods, precisely observed longitudes of all the equator cross- 
ings is sufficient for a harmonic determination from the simple drift regime equation. 

DETERMINATION OF H 22 AND H 33 FROM 24-HOUR SATELLITE DATA 

In the case of the 24-hour satellites, numerical integration has shown (Reference 11) that ad- 
ditional gravity perturbations to the regime equation (Equation 29a) are within the standard devia- 
tions of the nodal acceleration observations presented in Table 3. The record consists of 3 years 
of reduced equator crossing accelerations in 13 longitudinally separated free-drift arcs for Syncom 
2 and 3 and the Early Bird satellite, all of which have nearly circular orbits. The gravity coef- 
ficients C , S n in Table 3 are defined from the J , \ of Equation 11 as: 

nm 7 nm •'nm 7 nm * 


and 


s 


nm 


^nm 


sin m A. , 

nm 


to conform to another common way of expressing the potential as a fully linear series of harmonic 
functions. This representation, of less use physically, has the advantage of allowing a fully linear 
separation of harmonic effects convenient to the analysis and reduction of real data. The normal- 
ized longitude coefficients C 7 s are defined from c , s by 

0 nm 7 nm nm 7 nm J 



nm 



nm 
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Table 3 
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♦First Solution 
fSecond Solution 



The harmonics H 22 and H 33 reported in the right hand columns of Table 3 were derived by a 
weighted least-squares-solution of Equation 29a using the orbit data in the left columns of the table. 

The 24-hour satellite is in librational-resonance with longitude harmonics of all orders, the 
only restriction being that the difference between the order n and the longitude frequency m must be 
even. Nevertheless, the altitude of the satellite is so high that the effects of the third order har- 
monics are considerably damped with respect to that of the second order (h 22 ). This is the har- 
monic which reflects the earth's elliptical equator. 

The fairly extensive data in Table 3 was examined for sensitivity, first to the second order 
harmonic alone (results in parentheses in Table 3, and secondly to a combination of 2nd and 3rd 
order harmonics (results directly below those in parentheses). In the second solution, only small 
correlation was found between the effects of 2nd and 3rd order, and only minor changes from the 
2nd order coefficients of the first solution. This and additional analysis which include other relevant 
harmonics, appear to show that the sectorial harmonics of second and third order, h 22 and h 33 have 
been determined by these data alone to an absolute precision within the bounds shown in the second 
solution of Table 3. Comparison of these results with recent data from lower-altitude satellites 
shows reasonable agreement, except perhaps with respect to S 33 . 

The data tests for sensitivity to h 31 , h 42 and H 44 did not produce a significant change in the 
results for h 22 and h 33 or an improvement in their standard errors. It is concluded that the 24- 
hour data are not accurate enough yet or sufficiently widespread in longitude to allow a good de- 
termination of these other resonant harmonics through 4th order. 

DISCUSSION 

In mechanical systems, resonant phenomena occur when the periodicity of applied loads is 
close to the natural (unforced) periodicities of the system. Resonance in this sense is character- 
ized by greatly amplified vibrations of the system. Viewed internally, the applied loads are always 
in step, or in the same direction as the moment of the internal parts upon which they act, and so 
continually increase the vibrational energy of those parts. But viewed externally, the loads are 
purely sinusoidal, and so the mean translational energy of the center of mass of the system is un- 
changed. In the case of the satellite whose orbit yields a constant ground-track with respect to 
certain harmonics, the applied loads are not purely sinusoidal, but contain a bias component which 
drives the long-period libration. In the case of the constant ground-track orbit, resonance implies 
a change in the translational, or orbit energy of the satellite. The vibrational, or internal energy 
of the orbit can perhaps best be characterized by the eccentricity. In an eccentric orbit, the satel- 
lite can be imagined as vibrating continually about its mean anomaly and semimajor axis, or mean 
radius. This vibrational aspect of the orbit has been ignored by orbit-averaging the disturbing 
force. This is done because we are here interested only in long-term orbit energy changes. In 
terms of the mechanical resonance analogy, the full effect would be equivalent to a mass-spring 
assembly on wheels, vibrating and translating under the action of a biased, sinusoidally varying, 
external force. This paper has been concerned only with the steady, translational part of the 
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complete motion. In celestial mechanics, these changes are called secular, or effects which in- 
crease in time without apparent limit. In the case of commensurate orbits, the bias of the dis- 
turbing force will also change sinusoidally over a long period of time. This is the phenomenon 
called libration, having as its proper analogy the circular pendulum. (It should be noted again that 
the resonant orbits have been specified to be near circular. ) 

The "resonance" phenomenon that has been investigated by the orbit-averaging method here, 
might better be called secular or librational- circular pendulum effects. (In particular, we have 
really dealt only with the long term changes in the semimajor axis.) But the "in-step action" of 
the longitude harmonics on these commensurate orbits gives rise also to true vibrational reason- 
ances (in radius and anomaly), following the mechanical analogy (References 2, 18, and 19). 
Therefore, though the effect here analyzed is not commonly thought of as a resonance effect, it 
does arise in general for those orbits which also show true resonance effects; thus these orbits 
may be called resonant justifiably. When a "resonant gravity harmonic" is referred to in this 
paper, the word resonant is used only in the sense of being capable of producing the amplified 
changes of libration which are also properties of the resonant orbits. 

It has been pointed out that the true vibrational resonances in a commensurate orbit imply a 
buildup, or amplification, of the orbit eccentricity. This amplification would be without bound if 
true commensurability could be maintained. Of course it cannot, since the resonant orbit generally 
librates. In fact, as might be inferred from its dependence on the direction of the disturbing force, 
the eccentricity in a resonant orbit appears to go through the same long-period libration that the 
semimajor axis does. However, these true resonances of eccentricity (called "Dynamic Resonance" 
by Blitzer (Reference 19)) may not be as useful geodetically, since many more of the longitude har- 
monics H nm appear to contribute to them. 

The pendulum-like changes of the semimajor axis (describable from a separated solution of 
the coupled Equations 29 and 30) (References 13, 14), have a periodicity of years only for near- 
commensurate conditions. These are the librating orbits first discussed in Reference 20 for the 
24-hour case. In cases where the ground-track moves at a considerable rate, the basic circular 
pendulum Equations 29 and 30, with small modification (along the lines outlined in Reference 2), 
may still serve as the model of the evolution. In these cases of fast global circulation of the 
ground-track, evaluation of Equations 29 and 30 shows that the crossing rate and semimajor axis 
oscillate with a more rapid frequency but much diminished amplitude as compared to a libration 
case. For example, in a fast circulation case, the period is typically of the order of days, and the 
amplitude may be of the order of hundreds of meters (for the semimajor axis change). In full libra- 
tion cases, well away from the equilibrium positions, the amplitudes of the semimajor axis change 
may be of the order of tens of kilometers (Reference 13). The fast world-circulation regime (far 
from resonance) is actually a circulating pendulum phenomenon. The reduction in amplitude of the 
pendulum effect in a fast circulation regime is due to the rapidity with which the orbit-averaged 
bias force averages out over a global circulation. 

In spite of the small amplitudes, circulating pendulum resonance has been used to good effect 
in discriminating a number of very high order coefficients in the earth's field from a dense global 
tracking of the Navy's Doppler-Transit Satellites (Reference 4). 
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The limitations of the model used here to derive the key orbit evolution Equations 29 and 30 
should be re-emphasized. Though strictly speaking, these equations apply only to circular orbits of 
exact resonant period, they have been shown, by numerous examples (Reference 11) and in a rigorous 
way to be applicable with only slight correction or modification to a much wider class of satellites. 
This class includes those orbits with drift rates considerably exceeding the maximum permissible 
for libration as well as for noncircular orbits with eccentricities no higher than about 0.01. The 
equations, furthermore, give only the gross orbit-averaged effects due to the earth’s longitude 
gravity field. But over a period of the order of months, it may be presumed that, well away from 
any strong influence of the earth T s atmosphere or the moon T s gravity, these equations will describe 
the dominating long-term effect on the geographic configuration of all near-resonant, near -circular 
orbits, with the exception of inclination change due to the sun and moon. The extension of these regime 
equations to the class of eccentric resonant orbits is indicated by Gedeon, et. al. (See note at the bottom 
of page 15) 

SUMMARY AND CONCLUSIONS 

The long-term geographic evolution of circular-orbit satellite ground-tracks which are orig- 
inally near stationary has been found. The equation governing the motion of the near -stationary 
ground-track is essentially that of a circular pendulum for each relevant earth longitude harmonic. 
These relevant harmonics form a tenuous but infinite subset of the infinite set of longitude harmonics. 

For orbit periods very close to resonant, or rational with respect to the earth’s rotation rate, 
the sinusoidal ground-track librates with a maximum excursion of 360°/m, due to an individual 
relevant (resonant) earth longitude harmonic. The complete period of libration is two to ten years 
for the strongest resonances and depends strongly on the libration amplitude and inclination. 

The libratory-like evolution equations can be extended to far-from-stationary conditions to 
cover the smaller, but detectable, drift oscillations of higher frequency, called circulating pendulum 
resonance. 

Circular resonant orbits of periods 3, 4, 4.8, 6, 8, 12, 14.4, and 16 hours appear to be particu- 
larly suited to discriminating unambiguously the longitude harmonics H 98 , H 76 , H ss , H 54 , H 33 , H 32 , 

H 55 , and H 43 respectively, in a reasonably short period of time. As of 1966, the 24-hour resonant 
orbit satellites have provided the best discrimination of the two leading sectorial harmonics H 22 
and h 33 . 

The orbit determination for the resonant satellite used as a gravity probe would not need to be 
elaborate. It would have to insure that a reasonably circular orbit has been achieved. The inclina- 
tion must also be determined to about 0.1 degree accuracy. During free drift periods, the longitudes 
of all the equator crossings are sufficient data for a rapid harmonic determination. 

In conclusion, while the long-term librational resonances of the commensurate orbit satellites 
will be an annoying problem to the operators of communications satellites, they provide a powerful 
and unique means to solve many of the mysteries of the earth's complex gravity field. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, November 17, 1966 
188-48-01-05-51 
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Appendix 

Symbols 

The semimajor axis of the satellite's orbit (a* is the semimajor axis in units of 
earth radii). 

Phase and frequency constants of the disturbing function. 

Non-normalized cosine and sine gravity harmonic coefficients. 

Dominant, subdominant. 

The total energy (kinetic and two-body gravitational potential) of a satellite, 
e” , e = 2.718. . . 

A perturbation force (when superscripted with a prime a nondimensional force). 
When subscripted nmp, an inclination force function. 

An inclination force function, when superscripted with a prime, the real value of 
the force function. 

The tangential (along- track) component of the perturbation force F. 

Location of the Greenwich meridian. 

Signifying the gravitational harmonic term of order n and power m. 

Inclination of the satellite's orbit. 

The normalized amplitude, cosine coefficient and sine coefficient of the gravity 
harmonic H . 

nm 

The amplitude and phase of the non-normalized gravity harmonic H nm . 

A positive integer or zero. 

An inclination-dependent harmonic forcing function usually assumed constant 
over a libration period. 

As a subscript of in 8' (n, m, l) ; l = 1 or odd refers to a descending equator 
crossing, l = 2 or even refers to an ascending equator crossing. 
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The synodic period (in sidereal days) of the resonant orbit satellite, an integer. 
The number of integral days for the stationary ground track to begin to repeat 
itself. 

The associated Legendre polynomial of order n and degree m: 


p n m (x) 


(l -X 2 ) m/2 


d m P n (x) 

dx m 


where 


P (X S = ( 2n >- I x n __ n (n ~ 1) x n ~ 2 + n(n - 1) (n - 2) (n - 3) x n ~ 4 

n W 2 n (n!) 2 \ X (2 1 ) (1! ) (2n - 1) (2 2 ) (2! ) (2n - 1) (2n - 3) 

A positive integer, or zero. 

Number of orbits in the synodic period of the resonant satellite (n # sidereal days); 
an integer. 

The mean equatorial radius of the earth. 

The real part of { } . 

Position of the satellite (s 0 is its position at the start of the dynamics, and S' at 
the reference equator crossing). 

The period of the satellite’s orbit. 

Time from the beginning of the dynamics or some arbitrary zero. 

The minimum libration period of a resonant orbit. 

The gravitational potential function of the earth. 

The location of the vernal equinox. 

A small, but finite, change. 

Phase angles of the disturbing function. 

Argument of the ascending node, central angle in the orbit from the ascending 
node to the satellite. 

Generally, geographic longitude of the ascending equator crossing of the satellite. 
Can also refer to the geographic longitude of the descending equator crossing. 

An equilibrium longitude; stable, and unstable. 
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fi The earth T s Gaussian gravitational constant, 
cr Standard deviation (or estimate of it). 

4> Geographic latitude of the satellite. 

n The Keplerian element: Right Ascension of the Ascending Node, 
a) A frequency or angular rate. 

The earth’s inertial rotation rate. 
at Orbital revolution rate of the satellite. 

s 

(~ ) The r orbit average of the quantity ( ), except in reference to C~, S^, J nin . 

( ' ) d( )/dt , with the time increment in units of sidereal days. 

C *) d 2 ( )/dt 2 , with the time increment in units of sidereal days. 


NASA -Langley, 1967 30 


27 


